Identification and Expression Analysis of Putative Sugar Transporter Gene Family during Bulb Formation in Lilies.

Sugar transporters play important roles in plant growth and development, flowering and fruiting, as well as responses to adverse abiotic and biotic environmental conditions. Lilies (Lilium spp.) are some of the most representative ornamental bulbous flowers. Sugar metabolism is critical for bulb formation in lilies; therefore, clarifying the amount and expression pattern of sugar transporters is essential for further analyzing their roles in bulb formation. In this study, based on the transcriptome data of the Lilium Oriental hybrid 'Sorbonne' and Lilium × formolongi, a total of 69 and 41 sugar transporters were identified in 'Sorbonne' and Lilium × formolongi, respectively, by performing bioinformatics analysis. Through phylogenetic analysis, monosaccharide transporters (MSTs) can be divided into seven subfamilies, sucrose transporters (SUTs) can be divided into three subgroups, and sugars will eventually be exported transporters (SWEETs) can be divided into four clades. According to an analysis of conserved motifs, 20, 14, and 12 conserved motifs were predicted in MSTs, SUTs, and SWEETs, respectively. A conserved domain analysis showed that MSTs and SUTs contained a single domain, whereas most of the SWEETs harbored two MtN3/saliva domains, also known as a PQ-loop repeat. The LohINT1, which was predicted to have a smaller number of transmembrane structural domains, was cloned and analyzed for subcellular localization. It was found that the LohINT1 protein is mainly localized in the cell membrane. In addition, the expression analysis indicated that 22 LohMSTs, 1 LohSUTs, and 5 LohSWEETs were upregulated in 'Sorbonne' 1 day after scale detachment treatment, suggesting that they may regulate the initiation of the bulblet. A total of 10 LflMSTs, 1 LflSUTs, and 6 LflSWEETs were upregulated 4~6 months after sowing, which corresponds to the juvenile-to-adult transition phase of Lilium × formolongi, suggesting that they may also play a role in the accompanying bulb swelling process. Combined with quantitative real-time PCR (qRT-PCR) analysis, LohSTP8 and LohSTP12 were significantly overexpressed during the extremely early stage of bulblet initiation, and LflERD6.3 was significantly overexpressed during the growth of the underground bulblet, suggesting that they may be key sugar transporters in the formation of lily bulbs, which needs further functional verification.

LiNAC100 contributes to linalool biosynthesis by directly regulating LiLiS in Lilium 'Siberia'.

MAIN CONCLUSION: The transcription factor LiNAC100 has a novel function of regulating floral fragrance by directly regulating linalool synthase gene LiLiS. Lilium 'Siberia', an Oriental hybrid, is renowned as both a cut flower and garden plant, prized for its color and fragrance. The fragrance comprises volatile organic compounds (VOCs), primarily monoterpenes found in the plant. While the primary terpene synthases in Lilium 'Siberia' were identified, the transcriptional regulation of these terpene synthase (TPS) genes remains unclear. Thus, understanding the regulatory mechanisms of monoterpene biosynthesis is crucial for breeding flower fragrance, thereby improving ornamental and commercial values. In this study, we isolated a nuclear-localized LiNAC100 transcription factor from Lilium 'Siberia'. The virus-induced gene silencing (VIGS) of LiNAC100 was found to down-regulate the expression of linalool synthase gene (LiLiS) and significantly inhibit linalool synthesis. Conversely, transient overexpression of LiNAC100 produced opposite effects. Additionally, yeast one-hybrid and dual-luciferase assays confirmed that LiNAC100 directly activates LiLiS expression. Our findings reveal that LiNAC100 plays a key role in monoterpene biosynthesis in Lilium 'Siberia', promoting linalool synthesis through the activation of LiLiS expression. These results offer insights into the molecular mechanisms of terpene biosynthesis in Lilium 'Siberia' and open avenues for biotechnological enhancement of floral scent.

Transcriptome Analysis Reveals Coexpression Networks and Hub Genes Involved in Papillae Development in Lilium auratum.

Lilium is a genus of important ornamental plants with many colouring pattern variations. Lilium auratum is the parent of Oriental hybrid lilies. A typical feature of L. auratum is the presence of red-orange special raised spots named papillae on the interior tepals. Unlike the usual raised spots, the papillae are slightly rounded or connected into sheets and usually have hairy tips. To elucidate the potential genes regulating papillae development in L. auratum, we performed high-throughput sequencing of its tepals at different stages. Genes involved in the flavonoid biosynthesis pathway were significantly enriched during the colouration of the papillae, and CHS, F3H, F3'H, FLS, DFR, ANS, and UFGT were significantly upregulated. To identify the key genes involved in the papillae development of L. auratum, we performed weighted gene coexpression network analysis (WGCNA) and further analysed four modules. In total, 51, 24, 1, and 6 hub genes were identified in four WGCNA modules, MEbrown, MEyellow, MEpurple, and MEred, respectively. Then, the coexpression networks were constructed, and important genes involved in trichome development and coexpressed with anthocyanin biosynthesis genes, such as TT8, TTG1, and GEM, were identified. These results indicated that the papillae are essentially trichomes that accumulate anthocyanins. Finally, we randomly selected 12 hub genes for qRT-PCR analysis to verify the accuracy of our RNA-Seq analysis. Our results provide new insights into the papillae development in L. auratum flowers.

Lily (Lilium spp.) LhERF4 negatively affects anthocyanin biosynthesis by suppressing LhMYBSPLATTER transcription.

Anthocyanins are among the main pigments involved in the colouration of Asiatic hybrid lily (Lilium spp.). Ethylene, a plant ripening hormone, plays an important role in promoting plant maturation and anthocyanin biosynthesis. However, whether and how ethylene regulates anthocyanin biosynthesis in lily tepals have not been characterized. Using yeast one-hybrid screening, we previously identified an APETALA2 (AP2)/ETHYLENE RESPONSE FACTOR (ERF) named LhERF4 as a potential inhibitor of LhMYBSPLATTER-mediated negative regulation of anthocyanin biosynthesis in lily. Here, transcript and protein analysis of LhERF4, a transcriptional repressor, revealed that LhERF4 directly binds to the promoter of LhMYBSPLATTER. In addition, overexpression of LhERF4 in lily tepals negatively regulates the expression of key structural genes and the total anthocyanin content by suppressing the LhMYBSPLATTER gene. Moreover, the LhERF4 gene inhibits anthocyanin biosynthesis in response to ethylene, affecting anthocyanin accumulation and pigmentation in lily tepals. Collectively, our findings will advance and elucidate a novel regulatory network of anthocyanin biosynthesis in lily, and this research provides new insight into colouration regulation.

The antidepressive mechanism of Longya Lilium combined with Fluoxetine in mice with depression-like behaviors.

Traditional Chinese medicine is one of the most commonly used complementary and alternative medicine therapies for depression. Integrated Chinese-western therapies have been extensively applied in numerous diseases due to their superior efficiency in individual treatment. We used the meta-analysis, network pharmacology, and bioinformatics studies to identify the putative role of Longya Lilium combined with Fluoxetine in depression. Depression-like behaviors were mimicked in mice after exposure to the chronic unpredictable mild stress (CUMS). The underlying potential mechanism of this combination therapy was further explored based on in vitro and in vivo experiments to analyze the expression of COX-2, PGE2, and IL-22, activation of microglial cells, and neuron viability and apoptosis in the hippocampus. The antidepressant effect was noted for the combination of Longya Lilium with Fluoxetine in mice compared to a single treatment. COX-2 was mainly expressed in hippocampal CA1 areas. Longya Lilium combined with Fluoxetine reduced the expression of COX-2 and thus alleviated depression-like behavior and neuroinflammation in mice. A decrease of COX-2 curtailed BV-2 microglial cell activation, inflammation, and neuron apoptosis by blunting the PGE2/IL-22 axis. Therefore, a combination of Longya Lilium with Fluoxetine inactivates the COX-2/PGE2/IL-22 axis, consequently relieving the neuroinflammatory response and the resultant depression.

Lily LlHSFC2 coordinates with HSFAs to balance heat stress response and improve thermotolerance.

Heat stress transcription factors (HSFs) are core regulators of plant heat stress response. Much research has focused on class A and B HSFs, leaving those of class C relatively understudied. Here, we reported a lily (Lilium longiflorum) heat-inducible HSFC2 homology involved in thermotolerance. LlHSFC2 was located in the nucleus and cytoplasm and exhibited a repression ability by binding heat stress element. Overexpression of LlHSFC2 in Arabidopsis, tobacco (Nicotiana benthamiana), and lily, all increased the thermotolerance. Conversely, silencing of LlHSFC2 in lily reduced its thermotolerance. LlHSFC2 could interact with itself, or interact with LlHSFA1, LlHSFA2, LlHSFA3A, and LlHSFA3B of lily, AtHSFA1e and AtHSFA2 of Arabidopsis, and NbHSFA2 of tobacco. LlHSFC2 interacted with HSFAs to accelerate their transactivation ability and act as a transcriptional coactivator. Notably, compared with the separate LlHSFA3A overexpression, co-overexpression of LlHSFC2/LlHSFA3A further enhanced thermotolerance of transgenic plants. In addition, after suffering HS, the homologous interaction of LlHSFC2 was repressed, but its heterologous interaction with the heat-inducible HSFAs was promoted, enabling it to exert its co-activation effect for thermotolerance establishment and maintenance. Taken together, we identified that LlHSFC2 plays an active role in the general balance and maintenance of heat stress response by cooperating with HSFAs, and provided an important candidate for the enhanced thermotolerance breeding of crops and horticulture plants.

Cloning and Functional Characterization of 2-C-methyl-D-erythritol-4-phosphate cytidylyltransferase (LiMCT) Gene in Oriental Lily (Lilium 'Sorbonne').

2-C-methyl-D-erythritol-phosphate cytidylyltransferase (MCT) is a key enzyme in the MEP pathway of monoterpene synthesis, catalyzing the generation of 4- (5'-pyrophosphate cytidine)-2-C-methyl-D-erythritol from 2-C-methyl-D-erythritol-4-phosphate. We used homologous cloning strategy to clone gene, LiMCT, in the MEP pathway that may be involved in the regulation of floral fragrance synthesis in the Lilium oriental hybrid 'Sorbonne.' The full-length ORF sequence was 837 bp, encoding 278 amino acids. Bioinformatics analysis showed that the relative molecular weight of LiMCT protein is 68.56 kD and the isoelectric point (pI) is 5.12. The expression pattern of LiMCT gene was found to be consistent with the accumulation sites and emission patterns of floral fragrance monoterpenes in transcriptome data (unpublished). Subcellular localization indicated that the LiMCT protein is located in chloroplasts, which is consistent with the location of MEP pathway genes functioning in plastids to produce isoprene precursors. Overexpression of LiMCT in Arabidopsis thaliana affected the expression levels of MEP and MVA pathway genes, suggesting that overexpression of the LiMCT in A. thaliana affected the metabolic flow of C5 precursors of two different terpene synthesis pathways. The expression of the monoterpene synthase AtTPS14 was elevated nearly fourfold in transgenic A. thaliana compared with the control, and the levels of carotenoids and chlorophylls, the end products of the MEP pathway, were significantly increased in the leaves at full bloom, indicating that LiMCT plays an important role in regulating monoterpene synthesis and in the synthesis of other isoprene-like precursors in transgenic A. thaliana flowers. However, the specific mechanism of LiMCT in promoting the accumulation of isoprene products of the MEP pathway and the biosynthesis of floral monoterpene volatile components needs further investigation.

Regulatory mechanism of the miR172e-LbrAP2 module during the vegetative growth phase transition in Lilium.

MAIN CONCLUSION: It was proved for the first time that the miR172e-LbrAP2 module regulated the vegetative growth phase transition in Lilium, which provided a new approach to shorten the juvenile stage of Lilium, improved the reproduction rate, and reduced the propagation cost of Lilium commercial bulbs. Lilium is an ornamental bulb plant that takes at least 3 years to cultivate into commercial seed bulbs under natural conditions. The aim of this study was to shorten the Lilium expansion cycle. In this study, the growth cycle of lily tubers induced by low temperature of 15 °C was significantly shorter than that of tubers grown at a conventional temperature. Quantitative real-time PCR analysis showed that the expression patterns of miR172e and LbrAP2 were negatively correlated. GUS histochemical staining confirmed that miR172e and LbrAP2 in tobacco leaves interacted with each other after co-transformation. The shear sites of miR172e and its target gene, LbrAP2, upon binding, were identified by RLM 5' RACE analysis. In addition, miR172e and LbrAP2 showed opposite expression patterns after the transformation of Arabidopsis. miR172e overexpression accelerated the transition from juvenile to adult plants, whereas LbrAP2 overexpression inhibited this process, thus indicating that miR172e negatively regulated the target gene LbrAP2. Upregulation of the transcription factor LbrAP2 delayed the phase transition of plants, whereas miR172 inhibited the transcriptional translation of LbrAP2, thereby accelerating the phase transition. Low-temperature treatment of Lilium bulbs can shorten Lilium development, which provides a new approach to accelerating Lilium commercial bulb breeding and reducing breeding costs.

Low LdMYB12 expression contributes to petal spot deficiency in Lilium davidii var. unicolor.

Petal spots are widespread in plants, they are important for attracting pollinators and as economic traits in crop breeding. However, the genetic and developmental control of petal spots has seldom been investigated. To further clarify the development of petal spots formation, we performed comparative transcriptome analysis of Lilium davidii var. unicolor and Lilium davidii petals at the full-bloom stage. In comparison with the parental species L. davidii, petals of the lily variety L. davidii var. unicolor do not have the distinct anthocyanin spots. We show that among 7846 differentially expressed genes detected, LdMYB12 was identified as a candidate gene contributing to spot formation in lily petals. The expression level of LdMYB12 in the petals of L. davidii was higher than that in L. davidii var. unicolor petals. Moreover, overexpression of LdMYB12 led to the appearance of spots on the petals of L. davidii var. unicolor, accompanied by increased expression of anthocyanin synthesis-related genes. Taken together, these results indicate that abnormal expression of LdMYB12 contributes to petal spot deficiency in L. davidii var. unicolor.

Identification and functional analysis of LpNAC37 associated with somatic embryogenesis in Lilium pumilum DC. Fisch. based on transcriptome analysis.

Somatic embryogenesis (SE) is important for Lilium bulb propagation, germplasm conservation, and genetic transformation. The transition of somatic cells to embryonic cells is a critical step in SE, but the associated regulatory mechanisms have not been fully elucidated. Lilium pumilum DC. Fisch has a high regenerative capacity, and this study clarifies the critical timing of embryonic cell appearance in Lilium SE. Transcriptome sequencing using RNA-seq technology was performed on 5 representative samples from the early stage of Lilium SE. The 15 established cDNA libraries yielded 91.47 GB of valid data, and a total of 11,155 genes were consistently differentially expressed in the early stages of Lilium SE. GO annotation and KEGG pathway analysis of differentially expressed genes (DEGs) suggested that transcriptional regulation, hormone signaling, and stress response pathways play essential roles in the early stages of Lilium SE. WOX8, WOX11, SHR2, NAC37, AHP2, ANT, PIN1C, LAX2, LBD4, ACS12, YUC4, NFYB3, WRKY28, SAUR50, PYL9, and WRKY39 may be candidate genes for regulating early SE in Lilium. We further cloned LpNAC37, one of the key DEGs obtained from WGCNA and screening. LpNAC37 encodes a protein of 303 amino acids with a conserved NAM structural domain. The protein is a nuclear transcription factor with the highest homology to carrot DcNAC37. Overexpression of LpNAC37 suggested that LpNAC37 promotes embryonic callus formation in Arabidopsis. These results will help reveal the molecular mechanisms of the early stages of Lilium SE and advance the application of SE in Lilium propagation and genetic transformation.

A Modified Method for Transient Transformation via Pollen Magnetofection in Lilium Germplasm.

Lily (Lilium spp.) is a popular ornamental plant. Traditional genetic transformation methods have low efficiency in lily, thus development of a high-efficiency genetic transformation system is important. In this study, a novel transient transformation method involving pollen magnetofection was established and optimized pollen viability, and exogenous gene expression in magnetofected pollen and that of different germplasm were assessed. The highest germination percentage of Lilium regale pollen was 85.73% in medium containing 100 g/L sucrose, 61.5 mg/L H3BO3, and 91.5 mg/L CaCl2. A 1:4 ratio of nanomagnetic beads to DNA plasmid and transformation time of 0.5 h realized the highest transformation efficiency (88.32%). The GFP activity in transformed pollen averaged 69.66%, while that of the control pollen was 0.00%. In contrast to the control, transgenic seedlings obtained by pollination with magnetofected pollen showed strong positive GUS activity with 56.34% transformation efficiency. Among the lily germplasm tested, 'Sweet Surrender' and L. leucanthum had the highest transformation efficiency (85.80% and 54.47%), whereas L. davidii var. willmottiae was not successfully transformed. Transformation efficiency was positively correlated with pollen equatorial diameter and negatively correlated with polar axis/equatorial diameter ratio. The results suggest that pollen magnetofection-mediated transformation can be applied in Lilium but might have species or cultivar specificity.

LhANS-rr1, LhDFR, and LhMYB114 Regulate Anthocyanin Biosynthesis in Flower Buds of Lilium 'Siberia'.

The bulb formation of Lilium is affected by many physiological and biochemical phenomena, including flower bud differentiation, starch and sucrose accumulation, photoperiod, carbon fixation, plant hormone transduction, etc. The transcriptome analysis of flower buds of Lilium hybrid 'Siberia' at different maturity stages showed that floral bud formation is associated with the accumulation of anthocyanins. The results of HPLC-MS showed that cyanidin is the major anthocyanin found in Lilium 'Siberia'. Transcriptome KEGG enrichment analysis and qRT-PCR validation showed that two genes related to flavonoid biosynthesis (LhANS-rr1 and LhDFR) were significantly up-regulated. The functional analysis of differential genes revealed that LhMYB114 was directly related to anthocyanin accumulation among 19 MYB transcription factors. Furthermore, the qRT-PCR results suggested that their expression patterns were very similar at different developmental stages of the lily bulbs. Virus-induced gene silencing (VIGS) revealed that down-regulation of LhANS-rr1, LhDFR, and LhMYB114 could directly lead to a decrease in anthocyanin accumulation, turning the purple phenotype into a white color. Moreover, this is the first report to reveal that LhMYB114 can regulate anthocyanin accumulation at the mature stage of lily bulbs. The accumulation of anthocyanins is an important sign of lily maturity. Therefore, these findings have laid a solid theoretical foundation for further discussion on lily bulb development in the future.

Molecular and Metabolic Insights into Anthocyanin Biosynthesis for Spot Formation on Lilium leichtlinii var. maximowiczii Flower Petals.

Plants exhibit remarkable diversity in their petal colors through biosynthesis and the accumulation of various pigments. Lilium, an important cut and potted flower, has many coloring pattern variations, including bicolors and spots. To elucidate the mechanisms regulating spot formation in Lilium leichtlinii var. maximowiczii petals, we used multiple approaches to investigate the changes in petal carotenoids, spot anthocyanins, and gene expression dynamics. This included green petals without spots (D1-Pe and D1-Sp), yellow-green petals with purple spots (D2-Pe and D2-Sp), light-orange petals with dark-purple spots (D3-Pe and D3-Sp), and orange petals with dark-purple spots (D4-Pe and D4-Sp). D3-Pe and D4-Pe contained large amounts of capsanthin and capsorubin and small amounts of zeaxanthin and violaxanthin, which contributed to the orange color. In addition to cyanidin-3-O-glucoside, pelargonidin-3-O-rutinoside, cyanidin-3-O-rutinoside, and peonidin-3-O-rutinoside may also contribute to L. leichtlinii var. maximowiczii's petal spot colors. KEGs involved in flavonoid biosyntheses, such as CHS, DFR, and MYB12, were significantly upregulated in D2-Sp and D3-Sp, compared with D1-Sp, as well as in spots, compared with petals. Upregulated anthocyanin concentrations and biosynthesis-related genes promoted spot formation and color transition. Our results provide global insight into pigment accumulation and the regulatory mechanisms underlying spot formation during flower development in L. leichtlinii var. maximowiczii.

LiMYB108 is involved in floral monoterpene biosynthesis induced by light intensity in Lilium 'Siberia'.

KEY MESSAGE: We find that the MYB family transcription factor, LiMYB108, has a novel function to regulate the floral fragrance affected by light intensity. Floral fragrance determines the commercial value of flowers and is influenced by many environmental factors, especially light intensity. However, the mechanism by which light intensity affects the release of floral fragrance is unclear. Here, we isolated an R2R3-type MYB transcription factor LiMYB108, the expression of which was induced by light intensity and located in the nucleus. Light of 200 and 600 µmol m-1 s-1 significantly increased the expression of LiMYB108, which was consistent with the improving trend of monoterpene synthesis under light. Virus-induced gene silencing (VIGS) of LiMYB108 in Lilium not only significantly inhibited the synthesis of ocimene and linalool, but also decreased the expression of LoTPS1; however, transient overexpression of LiMYB108 exerted opposite effects. Furthermore, yeast one-hybrid assays, dual-luciferase assays, and electrophoretic mobility shift assays (EMSA) demonstrated that LiMYB108 directly activated the expression of LoTPS1 by binding to the MYB binding site (MBS) (CAGTTG). Our findings demonstrate that light intensity triggered the high expression of LiMYB108, and then LiMYB108 as a transcription factor to activate the expression of LoTPS1, thus promoting the synthesis of the ocimene and linalool, which are important components of floral fragrance. These results provide new insights into the effects of light intensity on floral fragrance synthesis.

Ethylene Response Factor LlERF110 Mediates Heat Stress Response via Regulation of LlHsfA3A Expression and Interaction with LlHsfA2 in Lilies (Lilium longiflorum).

Heat stress seriously affects the quality of cut lily flowers. The ethylene response factors (ERFs) participate in heat stress response in many plants. In this study, heat treatment increased the production of ethylene in lily leaves, and exogenous ethylene treatment enhanced the heat resistance of lilies. LlERF110, an important transcription factor in the ethylene signaling pathway, was found in the high-temperature transcriptome. The coding region of LlERF110 (969 bp) encodes 322 amino acids and LlERF110 contains an AP2/ERF typical domain belonging to the ERF subfamily group X. LlERF110 was induced by ethylene and was expressed constitutively in all tissues. LlERF110 is localized in the nucleus and has transactivation activity. Virus-induced gene silencing of LlERF110 in lilies reduced the basal thermotolerance phenotypes and significantly decreased the expression of genes involved in the HSF-HSP pathway, such as LlHsfA2, LlHsfA3A, and LlHsfA5, which may activate other heat stress response genes; and LlHsp17.6 and LlHsp22, which may protect proteins under heat stress. LlERF110 could directly bind to the promoter of LlHsfA3A and activate its expression according to the yeast one hybrid and dual-luciferase reporter assays. LlERF110 interacts with LlHsfA2 in the nucleus according to BiFC and the yeast two-hybrid assays. In conclusion, these results indicate that LlERF110 plays an important role in the basal thermotolerance of lilies via regulation of the HSF-HSP pathway, which could be the junction of the heat stress response pathway and the ethylene signaling pathway.

Lily HD-Zip I Transcription Factor LlHB16 Promotes Thermotolerance by Activating LlHSFA2 and LlMBF1c.

HD-Zip I transcription factors play important roles in plant development and response to abiotic stresses; however, their roles in thermotolerance are largely unknown. Through transcriptome analysis in lily (Lilium longiflorum), we isolated and identified a HD-Zip I gene differentially expressed at high temperatures, LlHB16, which belongs to the ß2 subgroup and positively regulates thermotolerance. The expression of LlHB16 was rapidly and continuously activated by heat stress. LlHB16 protein localized to the nucleus and exhibited transactivation activity in both plant and yeast cells, and its C-terminus contributed to its transcriptional activity. Overexpressing LlHB16 in Arabidopsis and lily improved thermotolerance and activated the expression of heat-related genes in both plants, especially that of HSFA2 and MBF1c. In addition, LlHB16 overexpression in Arabidopsis also caused growth defects, delayed flowering and abscisic acid (ABA) insensitivity. Further analysis revealed that LlHB16 directly binds to the promoters of LlHSFA2 and LlMBF1c and activates their expressions. Similarly, the expression of AtHSFA2 and AtMBF1c was also elevated in LlHB16 transgenic Arabidopsis lines. Together, our findings demonstrate that LlHB16 participates in the establishment of thermotolerance involved in activating LlHSFA2 and LlMBF1c, and LlHB16 overexpression resulted in ABA insensitivity in transgenic plants, suggesting that LlHB16 links the basal heat-responsive pathway and ABA signal to collaboratively regulate thermotolerance.

Expression Patterns of Key Genes in the Photoperiod and Vernalization Flowering Pathways in Lilium longiflorum with Different Bulb Sizes.

Lilium longiflorum is a wild Lilium, and its flowering transition requires a long period of cold exposure to meet the demand of vernalization. The responses of different sized bulbs to cold exposure and photoperiod are different, and the floral transition pathways of small and large bulbs are different. In this study, small and large bulbs were placed in cold storage for different weeks and then cultured at a constant ambient temperature of 25 °C under long day (LD) and short day (SD) conditions. Then, the flowering characteristics and expression patterns of key genes related to the vernalization and photoperiod pathways in different groups were calculated and analyzed. The results showed that the floral transition of Lilium longiflorum was influenced by both vernalization and photoperiod, that vernalization and LD conditions can significantly improve the flowering rate of Lilium longiflorum, and that the time from planting to visible flowering buds' appearance was decreased. The flowering time and rate of large bulbs were greatly influenced by cold exposure, and the vernalization pathway acted more actively at the floral transition stage. The floral transition of small bulbs was affected more by the photoperiod pathway. Moreover, it was speculated that cold exposure may promote greater sensitivity of the small bulbs to LD conditions. In addition, the expression of LlVRN1, LlFKF1, LlGI, LlCO5, LlCO7, LlCO16, LlFT1, LlFT3 and LlSOC1 was high during the process of floral transition, and LlCO13, LlCO14 and LlCO15 were highly expressed in the vegetative stage. The expression of LlCO13 and LlCO14 was different under different lighting conditions, and the flowering induction function of LlCO9 and LlFT3 was related to vernalization. Moreover, LlFKF1, LlGI, LlCO5, LlCO16, LlSOC1 and LlFT2 were involved in the entire growth process of plants, while LlCO6, LlCO16 and LlFT1 are involved in the differentiation and formation of small bulblets of plants after the inflorescence stage, and this process is also closely related to LD conditions. This study has great significance for understanding the molecular mechanisms of the vernalization and photoperiod flowering pathways of Lilium longiflorum.

WUSCHEL-related homeobox genes cooperate with cytokinin to promote bulbil formation in Lilium lancifolium.

The bulbil is an important vegetative reproductive organ in triploid tiger lily (Lilium lancifolium). Based on our previously obtained transcriptome data, we screened two WUSCHEL-related homeobox (WOX) genes closely related to bulbil formation, LlWOX9 and LlWOX11. However, the biological functions and regulatory mechanisms of LlWOX9 and LlWOX11 are unclear. In this study, we cloned the full-length coding sequences of LlWOX9 and LlWOX11. Transgenic Arabidopsis (Arabidopsis thaliana) showed increased branch numbers, and the overexpression of LlWOX9 and LlWOX11 in stem segments promoted bulbil formation, while the silencing of LlWOX9 and LlWOX11 inhibited bulbil formation, indicating that LlWOX9 and LlWOX11 are positive regulators of bulbil formation. Cytokinin type-B response regulators could bind to the promoters of LlWOX9 and LlWOX11 and promote their transcription. LlWOX11 could enhance cytokinin pathway signaling by inhibiting the transcription of type-A LlRR9. Our study enriches the understanding of the regulation of plant development by the WOX gene family and lays a foundation for further research on the molecular mechanism of bulbil formation in lily.

Comparative metabolome analysis unravels a close association between dormancy release and metabolic alteration induced by low temperature in lily bulbs.

KEY MESSAGE: The correlation between dormancy release and metabolic metabolic changes in lily bulbs during low temperature storage was investigated. Low temperature is a major environmental factor required for dormancy release in lily bulbs. Although great advances in plant metabolomics have been achieved, knowledge about the molecular basis of lily bulb metabolomes at different developmental stages in response to low temperature is still limited. In this work, the dormancy release, vegetative growth, flowering, metabolic profile and gene expression in the less dormant cultivar Lilium longiforum × Oriental hybrid 'Triumphator' (T) and the more dormant cultivar Lilium Asiatic hybrid 'Honesty' (H) were compared. Exposure to low temperature (LT) successfully promoted stem elongation, floral transition and flowering of both T and H bulbs. However, exposure to room temperature (RT) restricted stalk elongation of both T and H bulbs, and prohibited floral transition and flowering of H bulbs. Correspondingly, higher antioxidant enzyme activity and total primary metabolite contents were observed in the apical bud of T bulbs. Gene expression analysis revealed that expressions of LiFT, LiFLK, LiSOC1 and LiCBF were decreased, whereas the expression of LiSVP and LiFLC were increased, in the apical bud of H bulbs under RT storage condition. Our findings reveal that the growth and dormancy breaking of lily bulbs are closely associated with the metabolic changes in the apical buds during postharvest storage.

Overexpression of a novel heat-inducible ethylene-responsive factor gene LlERF110 from Lilium longiflorum decreases thermotolerance.

AP2/ERF (APETALA2/ethylene-responsive factor) family transcription factors are involved in various plant-specific processes, especially in plant development and response to abiotic stress. However, their roles in thermotolerance are still largely unknown. In the current study, we identified a heat-inducible ERF member LlERF110 from Lilium longiflorum that was rapidly induced by high temperature. Its protein was localized in the nucleus, and transcriptional activation activity was observed in yeast and plant cells. In addition, LlERF110 was able to bind to GCC- and CGG-elements, but not to DRE-elements. Overexpression of LlERF110 conferred delayed bolting and bushy phenotype, with decreased thermotolerance accompanied by a disrupted ROS (reactive oxygen species) homeostasis in transgenic plants. The accumulation of LlERF110 may activate certain repressors related to heat stress response (HSR) and indirectly damage the normal expression of heat stress (HS)-protective genes such as AtHSFA2, which consequently leads to reduced thermotolerance. Our results implied that LlERF110 might function as a heat-inducible gene but may hinder the establishment of thermotolerance.